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Abstract 

Forward-backward multiplicity correlations in symmetric collisions are cal- 
culated independently of the detailed form of the corresponding multiplicity 
distribution. Applications of these calculations to e^e~ annihilation and pp 
collisions confirm the existence of the weighted superposition mechanism of 
different classes of substructures or components. When applied to pp colli- 
sions in particular, clan concept and its particle leakage from one hemisphere 
to the opposite one become of fundamental importance. The increase with 
cm. energy of the correlation strength as well as the behaviour of the average 
number of backward particles vs. the number of forward particles are correctly 
reproduced. 
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1 Essentials on forward-backward multiplicity correlation in symmetric col- 
lisions 

The average number of charged particles generated in different events in the backward 
hemisphere (B), ub, is a function of the number of particles occurring in the forward 
hemisphere (F), rip, controlled by the correlation strength 6fb 

Cov[nB,nF] 

Ofb = ■ (1) 

y Var [n b] Var [n p] 

In hadron-hadron collisions |1|, ^ ^ Q the correlation strength parameter is rather large 
with respect to e~^e~ annihilation ^ and is growing with cm. energy in the total sample 
of events as shown in Table 0. 

In addition in e~^e~ annihilation at LEP energies it has been found |^ that 6fb ~ 
in the separate two- and three-jet sample of events. No information is available on the 
correlation strength in the separate samples of soft (no minijets) and semihard (with 
minijets) events in hadron-hadron collisions. 

2 The problem 

We want to calculate the parameter 6fb for the multiplicity distribution 

P{n)= Yl PtotAnp^nB), (2) 

nB+np=n 

where np and hb are random variables and Ptotai{nF,nB) is the joint distribution for the 
weighted superposition of different classes of events, i.e., 

^totai(^F, ns) = aPiinp, Ub) + (1 - a)P2{nF, Ub), (3) 
a being the weight of class 1 events with respect to the total. 



3 The general solution 

ab,Dl,{l + 62) + (1 - «)&2l?n,2(l + bi) + 

+ ia(l - a)(n2 - ni)2(l + 6i)(l + 62) 
aDl,il + b,) + il-a)Dl,il + b,) + 

+ ia(l - a){n2 - rii)2(l + bi){l + 62) 



Table 1: Experimental results on forward-backward correlation strength. 
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where bi are the correlation strengths of class 1 {i = 1) and class 2 {i = 2) events, Dn^i are 
the multiplicity distribution dispersions of class 1 (i = 1) and class 2 {i = 2) events and 
Hi the corresponding average charged multiplicity for class 1 {i = 1) and class 2 {i = 2) 
events. 

In case bi and 62 are zero (as in the separate two samples of events in e~^e~ annihilation) 
one finds 

^ -Q)(n2 -^1)^ 

aD.^ 1 + (1 - a)-D^,2 + - tt)(^2 - "-i)^ 

It should be pointed out that above formulas are independent from any specific form 
of the multiplicity distributions Pi and P2! They depend only on the weight alpha and 
average charged multiplicities and dispersions of the two classes of events. 



4 Applications of Eqs. (|^) and 

4-1 An intriguing application of Eq. to e+e^ annihilation 

Opal collaboration has found that forward backward multiplicity correlations are non 
existent in the separate two- and three-jet samples of events i.e. bi and 62 in the first 
general formula are zero and the correlation strength of the total sample of 2-jet and 3-jet 
events is equal to 0.103 ± 0.007. 

Using a fit to OPAL data with similar conditions to the jet finder algorithm for the 
separate samples of events we can determine all parameters in formula @ and test its 
prediction with the experimental finding. 

It turns out that the values of the parameters]^ needed in are a = 0.463, ni = 18.4, 
77,2 = 24.0, Df = 25.6, D2 = 44.6 and the predicted value of fepB is 0.101, in extraordinary 
agreement with experimental data! 



4-2 A suggestive application of Eq. to pp collisions 

The application of (|]) to pp collisions leads to unsatisfactory results but opens a new per- 
spective: forward-backward multiplicity correlations cannot be neglected in the separate 
components. Accordingly Equation and not (|]) should be used. 

Repeating the same approach done in e'^e~ annihilation for calculating bp-Q (i.e., as- 
suming that in the separate samples of events FB multiplicity correlations are absent, 
61 = 62 = 0) in the case of pp collisions at 546 GeV cm. energy and using Fuglesang's 
fit|^ to soft and semihard events (accordingly a = 0.75, ni = 24.0, 7x2 = 47.6, D"^ ^ = 106, 

D^2 = 209) one finds 6fb = 0.28 [b^^^ = 0.58). The theoretical prediction in this case 
is too small! It is clear that our working hypothesis was not correct in this case. In 
conclusion forward-backward multiplicity correlations are needed in each class of events, 
i.e., bi and 62 should be different from zero, and after their determination general formula 
(|) and not formula (^ should be used! 

Results in 4.1 and 4.2 are a striking test of the existence of the weighted superposition 
effect, only a guess up to now. 
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5 A new theoretical problem 



Following above conclusions the next problem is how to determine bi and 62 when explicit 
data on forward-backward multiplicity correlations in the two separate samples of events 
are lacking and 6fb of the total sample is known from experiments. 

The generality of Equation (^) should be limited by introducing additive assump- 
tions inspired by our phenomenological knowledge of the particle emission process in the 
collision under examination. 

Assuming for instance that 

a. particles are independently produced in the collision, 

b. binomially distributed in the forward and backward hemispheres, 
it is found that 

D"^ ■ - n- 

' ' (6) 



where i5„ j and are the dispersion and the average charged multiplicity of the overall 
multiplicity distribution of each component being as usual i = 1,2. 
Assuming next that 

c. the multiplicity distribution in each z-component is NB (Pascal) with parameters 
and ki (an assumption which is suggested by the success of the weighted superposition 
mechanism of NB(Pascal)MD's in describing shoulder effect in charged particle multiplic- 
ity distributions and Hq vs q oscillations and which we hardly would like to abandon), 

we find 

= iirk- 

Accordingly bi values can be calculated by using again Fuglesang's fit parameters on 
the two components at 546 GeV cm. energy. After inserting in the general formula (^) 
these parameters we find 6fb = 0.78. 

A too large value with respect to the experimental one (6fb = 0.58)! This result leads 
to the following question: Which one of above mentioned apparently quite reasonable 
assumptions should be modified? 

Our guess is that charged particle FB multiplicity correlation is not compatible with 
independent particle emission but is compatible with the production in cluster, i.e., clan 
within a NB (Pascal) MD framework. An idea which we propose to develop and to explore 
in the following. 



6 Clan concept is of fundamental importance 

Successive steps of our argument ||T0| are 

i) the joint distribution -Ptotai(^F, 'T'b) is written as the convolution over the number 
of produced clans and over the partitions of forward and backward produced particles 
among clans: 

P(n^,nB) = J] P(Ap,Ab) Y1 PF{m'p,m'B\NF)pB{m'^,m'i,\NB). (8) 

Np,Nb m'p+m'p=np 
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ii) forward backward hemispheres symmetry property is used 



PF{n,m\N) = pBirrijnlN). (9) 

iii) leakage parameter p is introduced: it controls the probability that a binomially 
distributed particle generated by one clan lying in one hemisphere has to leak in the 
opposite hemisphere, q is the leakage parameter working in the symmetric direction, 
p + q = 1 (notice that p = 1 or g = means no leakage, the variation domain of p is 
0.5 < p < 1 and when p < 0.5 the clan is classified in the wrong domain). 

iv) covariance 7 = {{fip — P'F){f^B — fis)) oi hf forward and hb backward particles 
within a clan for forward and backward binomially distributed particles generated by 
clans is also introduced. 

v) clans are binomially produced in the forward and backward hemispheres with the 
same probability and particles within a clan are independently distributed in the two 
hemispheres. 

It follows for each z-component 

Dl- A{d%^{N)){p - qf + AN^/nl 
D% + A{dl^{N)){p - qf - AN^/nl + 2NDl/nl 
^ Dl/n - Dl/n, - A{dl^{N)){p - qfnjN + 47/n, 
Dl/n + Dl/n, + 4«^(iV))(p - qfnjN - A^/n,- 

Eq. (0) assuming NB (Pascal) behavior with characteristic Ui and fcj parameters for 
each component, binomial clan distribution in the two hemispheres, binomial distribution 
in the two hemispheres of logarithmically produced particles from each clan according to 
clan structure analysis gives 

2niPiqi ^^^^ 



rii + ki- 2niPiqi 



Accordingly the problem is therefore reduced to determine leakage parameters Pi in 
the two classes of events! 

Notice that in the limit nj — > 00, for decreasing fcj, bi depends on pi only. 



7 A phenomenological argument for determining leakage parameters pi 

By assuming that the semihard component is negligible at 63 GeV cm. energy and know- 
ing b-pB from experiment at such energy, equation ([Tl|) allows to determine psoft (0.78); the 
relatively small variation of nc,soft from 63 GeV to 900 GeV (it goes from 2 to ~ 2.44) 
leads to the conclusion that the leakage parameter for the soft component Psoft can be 
considered in the GeV domain nearly constant, i.e., Psoft = 0.78: therefore the correlation 
strength for the soft component at 546 GeV cm. energy, 6soft(546 GeV), can easily be 
determined. 

The germane equation for 6semihard(546 GeV) contains of course the unknown param- 
eter ^semihard at the cm. energy of 546 GeV. By inserting in equation (^) for fopB (total) 
^soft(546 GeV) = 0.78 and &semihard(546 GeV) as given by equation (0) with unknown 



5 



1 

b 



1 1 1 1 M M I 

soft 


1 1 1 M 1 I 1 1 1 M 1 M I . ■ 


_ semihard 




- total 






/ — — ""(D 
^ — — 

/ 


1 1 1 1 1 1 1 1 1 


/ 

/ 

/ 

^ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



10 10^ 10^ 10"^ 

cm. energy (GeV) 



Figure 1: Predictions for tlie correlation coefRcients for each component (soft and semi- hard) and for 
the total distribution in collisions in scenario 2. Three cases are illustrated, corresponding to the three 
numbered branches: leakage increasing with -Js (upper branch, ®), constant leakage (middle branch, @) 
and leakage decreasing with -y/s (lower branch, (D). Leakage for the soft component is assumed constant 
at all energies. The dotted line is a fit to experimental values. 



^semihard parameter, Psemihard at 546 GeV can be calculated from the experimental value 
of 6fb (total) = 0.58. It is found Psemihard(546 GeV) = 0.77. 

Since r?,c,semihard does not vary too much in the GeV region (it goes from 1.64 at 200 
GeV cm. energy to 2.63 at 900 GeV cm. energy, a relatively small variation which will 
hardly affect the corresponding leakage parameter in this domain) it is not hazardous to 
take ^semihard ~ coustaut in the same region. 

Under just mentioned assumptions 

a. the correlation strength cm. energy dependence is correctly reproduced in the GeV 
energy range from ISR up to UA5 top cm. energy and follows the phenomenological 
formula 6fb = -0.019 + 0.061 Ins (see Fig. |I|). 

h. when extrapolated to the TeV energy domain in the scenarios discussed in Ref. 
with the same values of psoft obtained in the GeV region (nc,soft(14 TeV) being 2.98 
makes this guess acceptable) and ^semihard also constant (a too strong assumption of 
course), a clean bending effect in fopB vs. Ins is predicted. Bending effect is enhanced 
or reduced by allowing Psemihard to increase (less leakage from clans and more bending) or 
to decrease logarithmically with cm. energy (more leakage from clans and less bending). 
Energy dependence of leakage parameter for the semihard component is clearly expected 
in the TeV region in a scenario with strong KNO scaling violation in view of the quite 
large average number of particles per clan with respect to that found at 900 GeV cm. 
energy (it goes from 2.63 at 900 GeV up to 7.36 at 14 TeV). See again Fig. ||. 

c. in addition nsinp) behavior at 63 GeV cm. energy (ISR data) is quite well de- 
scribed in terms of the soft component (single NB) only and at 900 GeV cm. energy 
(UA5 data) in terms of the weighted superposition of soft and semihard components, i.e., 
of the superposition of two NB(Pascal)MD's. (See Fig. ||, where the second case is shown). 
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(b) 900 GeV |ri|<4 




Figure 2: Results of our model for UBinp) vs. np compared to experimental data in the pseudo-rapidity 

interval < 4 at 900 GeV. 



8 Conclusions 

Weighted superposition mechanism of two samples of events describes forward backward 
multiplicity correlations in e^e' annihilation independently of the specific form of the 
charged particle MD in the different classes of events: only the average numbers of particles 
and related dispersions in addition to the weight factor are needed. 

In order to describe forward backward multiplicity correlations in pp collisions lack of 
information on FB multiplicity correlations in the separate components is demanding to 
specify the form of particle multiplicity distributions of the two components. 

The choice of NB(P)MD for each component (supported by its success in describ- 
ing shoulder effect and Hq vs q oscillations) outlines the role of clan properties in this 
framework and allows to determine correctly fepB energy dependence for the total sample 
of events in the GeV region. Its bending in the TeV region within possible scenarios 
discussed in the literature is predicted. 

nBinp) vs Up trend is also nicely reproduced at 63 GeV (only soft component is 
assumed to contribute) and at 900 GeV (superposition of soft and semihard components 
is used), and its behavior in the TeV energy range predicted. 

Last but not least we have found that our study on FB multiplicity correlations in 
pp collisions when extended to the TeV energy region assuming KNO scaling violation 
for the semihard component enhances the intriguing connection already shyly anticipated 
in the GeV region between particle populations within clans, particle leakage from clans 
in one hemisphere to the opposite one and the superposition effect between different 
components. Clan concept appears in this framework as a powerful tool which goes far 
beyond its simple statistical interpretation and raises the question on its real physical 
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significance: an interesting but also compulsory question for future experimental work in 
pp collisions and not only. 
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